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Abstract 

The question whether a spin-imbalanced Fermi gas can accommodate the Fulde-Ferreh-Larkin- 
__i _ Ovchinnikov (FFLO) state has been the subject of intense study. This state, in which Cooper 

^ ■ pairs obtain a nonzero momentum, has hitherto eluded experimental observation. Recently, we 



demonstrated that the FFLO state can be stabilized in a 3D Fermi gas, by adding a ID periodic 
potential. Until now it was assumed that the FFLO wave vector always lies parallel to this periodic 



'^ ■ potential (FFLO-P). In this contribution we show that, surprisingly, the FFLO wave vector can 

O ■ also lie skewed with respect to the potential (FFLO-S). Starting from the partition sum, the saddle- 

point free energy of the system is derived within the path-integral formalism. Minimizing this free 

(N 

^ ■ energy allows us to study the different competing ground states of the system. To qualitatively 

m" 

understand the underlying pairing mechanism, we visualize the Fermi surfaces of the spin up and 

^N . spin down particles. From this visualization, we find that tilting the FFLO wave vector with respect 

»■. 

i^ \ to the direction of the periodic potential, can result in a larger overlap between the pairing bands 

of both spin species. This skewed FFLO state can provide an additional experimental signature 
, , for observing FFLO superfluidity in a 3D Fermi gas. 
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I. INTRODUCTION 



Over the last decade, immense progress has been made in the field of the physics of ultra- 
cold atoms. These systems provide a versatile tool for studying various quantum many-body 
phenomena [l|, l2[]. One example of the major breakthroughs in this field is the realization 
of a molecular Bose- Einstein condensate (BEG) (3|, \^ and the observation of pair formation 
„ear a Feshbach reso„a.ce ffl in a strongly interacting Fermi gas. These experinrenta. 
discoveries preluded the observation of superfiuidity in a gas of ultracold fermions [lOj. The 
advantage of using ultracold Fermi gases to study superfiuidity is that they allow to explore 
a wide range of parameter space. This is due to the tunability of the system parameters, 
such as the interaction strength and the spin imbalance. The possibility of adapting the 
interaction strength through the use of a Feshbach resonance ll| has allowed to study su- 
perflu,d pairing in the crossover from a Bardeen-Cooper-Sdrr.effer (BCS) state of weakly 
bound Cooper pairs to a BEG of tightly bound molecules p^| . Furthermore, in a Fermi gas 
composed of a mixture of two hyperfine states (labeled "spin up" and "spin down"), the 
ratio between the number of atoms in these two different states can be controlled with great 
precision. This last achievement has provided a unique experimental tool to study the effect 
of population imbalance on superfiuidity. The first theoretical contribution in this context 
was provided by Glogston and Ghandrasekhar, who predicted that a first order transition 
:'rom a superfiuid to a normal interacting Fermi gas would occur at a critical polarization 
13| . This quantum phase transition was indeed observed experimentally IJ, ll5|, along 



with the fact that when the system is in the BGS superfiuid state, the excess particles are 
expelled from this state, leading to phase separation. 

The question that has emerged here is whether there exist other, exotic quantum many 
body states that allow polarized superfiuidity. In 1964, a new state was introduced inde- 
pendently by Fulde and Ferrell 16| and by Larkin and Ovchinnikov 17|]. They proposed 
that a spin-polarized superfiuid can be formed by creating Gooper pairs with finite center- 
of-mass-momentum. In the last seven years, extensive theoretical research has been done 
on this exotic superfiuid state [18|. In ID and quasi- I D sy stems, the Fulde- Ferrell-Larkin- 



Ovchinnikov (FFLO) state is predicted to be stable 19|, and recently the first indirect 
experimental evidence for FFLO in a ID system has been found ^0|. In general, lower 
dimensionality favors the FFLO state because of Fermi surface nesting at the FFLO wave 



vector. For the case of a three dimensional (3D) Fermi gas, it was found that the FFLO 



state only occurs on a restricted area of the BCS-to-BEC-crossover phase diagram 2l|, |22 |. 



Up till know, the FFLO state has eluded experimental observation in a 3D Fermi gas. As 



an attempt to pave the path towards its experimental disco very , t 



le idea was proposed to 



enhance FFLO pairing by using a 3D periodic optical lattice J23|, l2J|. Recently, we proposed 
to subject the 3D Fermi gas to a one dimensional (ID) periodic potential as an alternative 
to achieve this goal [25] . The ID potential introduces an asymmetry into the system, which 
results in an energetically preferred direction for the FFLO wave vector, thus lowering the 
energy of this state compared to the normal state and to the BCS state. 

In our previous work, it was assumed that the FFLO wave vector lies parallel to the 
periodic potential, since this is the energetically preferred direction. In this paper we show 
that, contrary to this intuitive expectation, the wave vector of the FFLO state can also lie 
skewed with respect to the direction along which the periodic potential lies. To qualitatively 
understand this counterintuitive phenomenon, we present a visualization of the underlying 
pairing mechanism by plotting the Fermi surfaces of the spin up and spin down fermions. 
This will pinpoint the effect of the periodic potential on the FFLO pairing mechanism. Our 



paper is structured as follows. In section [IT 
the use of the path integral formalism 



m 



we calculate the free energy of the system through 



27j . By choosing an appropriate saddle point, we 



incorporate the possibility of the FFLO state in our description, where we allow the FFLO 
wave vector to lie in an arbitrary direction. In section UlIAl we minimize the free energy and 
discuss the different competing ground states of the system. Subsequently, in section IIIIBI 
we construct the phase diagram as a function of the total and imbalance chemical potential 
and in section IIIICI we visualize and discuss the changes in the pairing mechanism of the 
FFLO state due to the presence of the periodic potential. Finally in section IIVI we draw 
conclusions. 



II. PATH INTEGRAL TREATMENT 

The starting point of our analytic treatment is the partition sum of a 3D spin-imbalanced 
Fermi gas, written in path integral form in units h = 2m = 1 

Z = VtpVip exp I - ^ ^ ^k,u;„,a [-iuJn + e(k±,k^) - jj,^] ^k,, 

"^ \ k,n cr 
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k,n q,m / 

with /3 = l/ksT the inverse temperature and V the volume of the system. In expression 
([1]), the fermionic fields are described by two Grassmann variables ip and ip- Furthermore, 
k = (k_|_, A;^) is the single-particle momentum and Un = (2n + l)7r//3 and Qm = ^rmr/P 
represent the fermionic and bosonic Matsubara frequencies respectively. In the single-particle 
term, e (k^, k^) is the energy dispersion and the chemical potential fi^ fixes the number of 
spin up (a =f) and spin down {a =\.) particles. In ([T]), the interaction between fermions is 
modeled using a delta-function pseudo potential, given by V (r, r') = g6 {r ~ r'), where g is 
the renormalized interaction strength, which is given by 

1 1 r dk 1 



Svr J {2n)^26{k^,k,y ^'^^ 

where we have set the unit of length to a^ = — 1. We set the scattering length to a negative 
value because our aim is to study the FFLO state, which, to the best of our knowledge, 
only occurs on the BCS side of the BCS-to-BEC crossover. In this part of the crossover, the 
scattering length is always negative, hence the choice of a^. 

Our goal is to study the pairing mechanism of the FFLO state under the influence of a 
ID periodic potential. Imposing such a potential on the system results in a change in the 
energy dispersion. Here we assume that the potential is deep enough so that the dispersion 
is of the following tight-binding form 

irk 



e{k^,k,) = ki + S 



1 — cos 



(3) 



Ql 

where Ql = 2tt/X is the wave vector of the periodic potential with A the wavelength of 
the periodic potential, and 6 is the bandwidth. In expression ([3]), the periodic potential is 
assumed to lie along the z-direction. This convention will be utilized throughout the rest of 
this paper. Physically, the tight-binding approximation corresponds to the situation where 



the overlap between wave functions of particles on neighboring sites is taken into account 
and the next-to-nearest neighbor hopping is neglected. For a potential depth Vq > 4:Er, 
he tight-binding dispersion (^ agrees to within less than five percent with the exact result 



281], where Eji is the recoil energy which in our units is given by E^i = (27r/A)^. In ([3]) the 
bandwidth 6 is a. function of Eji and of the potential depth Vq. An analytic expression for S 
can be derived by solving the ID Mathieu equation in the limit Vq >> Er [29]. The result 



is given by 




^ = 4 ('^Vexp 1-2^/^1. (4) 



To calculate the partition sum ([T]) we use the Hubbard-Stratonovich transformation, 
which introduces two auxiliary complex bosonic fields A and A and simultaneously re- 
duces the fourth order interaction term to two second order interaction terms: ip^ipiilJiilj^ — >■ 
Aip^ipl+ilj^ilj^A. This transformation takes into account solely the Bogoliubov channel and is 
exact as long as there are no other competing channels that have comparable contributions. 
For the purpose of this paper, which focuses on the BCS side of the BCS-BEC crossover and 
on zero temperature, the use of Hubbard-Stratonovich is justified. Recently, an alternative 



approach to Hubbard-Stratonovich has been introduced 30|] based on Feynman-Kleinert 



variational perturbation theory 3l|. This treatment, however, lies beyond the scope of the 
present paper. 

Subsequently, the saddle-point approximation is introduced, by which the bosonic path 
integral is reduced to the single most contributing term 

JvA exp [S (A)] ^ exp {-S [A,, (q)]} , (5) 

where Agp, which in general depends on q, is chosen so that it minimizes the action S. 
To describe the FFLO state, we take a specific form for Asp so that the bosonic pairs are 
allowed to have a finite center-of-mass momentum Q 



A.p(q) = V/5^(5q,QA, (6) 

with 

Q = (Q±,Q.). (7) 

where the factor \^PV ensures that A has units of energy. The two variational parameters 
A and Q are interpreted respectively as the band gap of the system (or, equivalently, as 



the binding energy of the bosonic pairs) and the wave vector of the FFLO state. Here, we 
allow Q to have a nonzero perpendicular component Q^. Since the system exhibits an axial 
symmetry around the z-axis, it suffices to determine the magnitude |Q^| of this perpendic- 
ular component. After applying the saddle-point approximation, the only remaining path 
integral is Gaussian and can be calculated exactly. Finally, the Matsubara summation can 
be performed analytically, which results in the following expression for the saddle-point (sp) 
free energy Qgp that reads, in the zero temperature limit 

r.27r 
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dk\ k\ 
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where the following shorthand notations were used 



Ck,Q = C - 5 sin 

-E'k,Q = 



1 — COS 
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Here we have introduced /i = (/i-|- + /i^) /2 and ( = (/i^ — /i^) /2 which represent the total and 
imbalance chemical potential, respectively. To the best of our knowledge it is impossible to 
calculate expression (|8]) completely analytically. However, this expression can be simplified 
to a certain degree by integrating out the radial component k±. The details of this calculation 
are given in appendix |X1 

III. RESULTS AND DISCUSSION 

A. The different competing ground states 



To determine the ground state of the system, we minimize the saddle point free energy 
dH]) with respect to the three variational parameters (A,^^, |Q±|), for given values of the 
two thermodynamic variables /i and (. The values of these three variational parameters in a 
given minimum determine the ground state of the system, which can be any of the following 



four states: 

^ BCS ^ A ^ 0, Q, = 0, IQxl = 

FFLO-P ^ A ^ 0, Q, 7^ 0, IQxl = 

FFLO-S ^ A 7^ 0, Q, 7^ 0, IQxl ^ 

Normal — )• A = 

where BCS signifies the spin-balanced superfluid state. Here we make a distinction between 
the FFLO-P and the FFLO-S state, which respectively denote that the FFLO wave vector 
lies parallel (|Q_l| = 0) or skewed (|Q_l| 7^ 0) with respect to the direction of the periodic 
potential. The main question addressed in this paper is whether the FFLO-S state can be 
the ground state of the system, given certain values of fi and (. To show the competition 
between the different ground states it is instructive to look at a contour plot of the free 
energy as a function of A and Qz, as shown in Fig. [T] The third parameter |Q_l| is held 
constant in a given plot. In Fig. [H the various plots were made at the same value of n 
{fi = 4.44) and for increasing values of (. In Fig. [1] (a), for the lowest value of (, the system 
is in the spin-balanced BCS state. This state is characterized by a nonzero band gap A, 
and by zero momentum Q of the fermionic pairs. Figure [1] (a) shows that at finite values 
of (, the Fermi gas can still be spin-balanced. This is because the nonzero binding energy 
of the bosonic pairs has to be overcome in order to break up these pairs. The imbalance 
chemical potential ( can be interpreted as a Zeeman energy which tries to align the spins 
in a preferred direction. When ( increases, this magnetic energy increases likewise and will 
eventually become larger than the binding energy of the Cooper pairs, at which point the 
system will make a transition into a polarized state. This transition can be seen in Fig. 
[H (b), where a second minimum emerges that starts to compete with the BCS minimum. 
This new minimum lies at (A 7^ 0, Qz 7^ 0) but still at |Q_l| = and is hence the FFLO-P 
state, which has a wave vector parallel to the periodic potential. The transition from BCS 
to FFLO-P is of first order, because there is a jump in the value of the band gap A. When ( 
increases further, one in principle expects the FFLO-P state to continuously go over into the 
normal state, as the band gap A will gradually decrease to zero when ( increases. We indeed 
find this behavior at low values of /i (/i ^ 3.0). However, in the case of Fig. [1] (at yU = 4.44) 
the system behaves differently. When taking a closer look at the FFLO-P minimum in Fig. 
[1] (b), we see that the normal minimum (at A = 0) starts to compete with the FFLO-P 
minimum, which seems to imply a first order (FFLO-P)-to-normal transition (see Fig. [1] 
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FIG. 1: Several contour plots of the free energy as a function of the band gap A and the z- 
component of the FFLO wave vector Q^- The value of the perpendicular component |Q^| is held 
constant in each plot. Darker regions represent lower values of the free energy. Four different ground 
states can be identified: (a) a spin-balanced BCS superfluid (A 7^ 0, Qz = 0, |Q_l| = 0), (b)/(c) the 
FFLO-P state (A / 0, Q^ / 0, |Q^| = 0) , (e) the FFLO-S state (A / 0,Q^ / 0, |Q_l| / 0), and (f) 
the normal state (A = 0). In figures (c) and (d), the emergence of a local normal minimum that 
competes with the FFLO-P minimum is shown. There will however be no first order transition 
from FFLO-P to normal (c) — )■ (d), because at the transition value of ^ (C = 1.0) where the normal 
minimum lies lower than the FFLO-P minimum (d), a new global minimum emerges at |Q^| 7^ 
(e). This means that the system becomes an FFLO-S state. The behavior shown in figures (c), 
(d) and (e) only occurs for ^ ^ 3.0. For lower values of /U the system makes a transition from BCS 
over FFLO-P to Normal. In all contour plots we used /u = 4.44, Vq = 6, and A = 1200 nm. 

(c) and (d)). This is not what really happens though, because at the transition value of C, 
((^ = 1.0) where the normal minimum lies lower than the FFLO-P minimum (Fig. [11(d)), a 
new global minimum emerges at |Q_l| 7^ (Fig. [11(e)), meaning that the system has made 
a transition from the FFLO-P state into the FFLO-S state. This shows that, at sufficiently 
high values of the total chemical potential /i, the system can favor an FFLO state with 
a wave vector that lies skewed with respect to the direction of the ID periodic potential. 
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When ( is increased further, the value of A for the FFLO-S state goes continuously to zero, 
and the system makes a transition into a normal interacting Fermi gas (Fig. [11(f)). 

B. The phase diagram as a function of the chemical potentials 

The examples shown in Fig. [1] were only for one particular value of the total chemical 
potential fi. When the free energy is minimized and therefore the ground state is determined 
for a large set of values for /i and (, we obtain the phase diagram of the 3D imbalanced 
Fermi gas in the presence of a ID periodic potential at zero temperature, as shown in Fig. 
[2j This figure shows that the FFLO-S state can indeed be the ground state of the system. 
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FIG. 2: Phase diagram of a spin-imbalanced 3D Fermi gas in the presence of a ID periodic potential, 
as a function of the total and imbalance chemical potentials /x and (^. The skewed FFLO state 
(FFLO-S) can be the ground state of the system but only at higher densities, when the Fermi 
energy lies above the energy corresponding to the top of the lowest Bloch band: Ep > 25. The 
region of the FFLO-P state (wave vector parallel to the periodic potential) is divided into two 
parts, where the z-component of the wave vector of FFLO Qz is either smaller or equal to the laser 
wave vector Qp. The total FFLO region (FFLO-P + FFLO-S) is significantly larger than in the 
known 3D case without periodic potential. In this figure we used Vq = 6 and A = 1200 nm for the 
depth and wavelength of the periodic potential respectively. 



Furthermore, this state only occurs at higher densities, more specifically when the Fermi 
energy lies above the energy corresponding to the top of the first Bloch band: Ep > 2S. 
The FFLO-P region can be subdivided into two regions. In one region, the z-component 
Qz of the wave vector of the FFLO-P state is smaller than the wave vector of the periodic 
potential Ql. In the other region, this z-component equals the wave vector of the potential. 
The FFLO-S state only occurs when Qz = Ql- This observation along with the fact that 
the FFLO-S state only occurs when Ep > 25 will be explained in the next subsection, where 
we illustrate the FFLO pairing mechanism in momentum space. As we have shown before 



25|, a striking feature of the phase diagram in Fig. |2]is the significant increase in the area of 
the FFLO state, compared to the case of a 3D Fermi gas without periodic potential. In the 
latter case, the FFLO area would barely be visible on the scale of Fig. |2i This is a result of 
the fact that the periodic potential lowers the energy for the formation of the FFLO state, 
relative to the BCS and the normal state. 

C. The influence of the periodic potential on the pairing mechanism of the FFLO 
state 

The occurrence of the FFLO-S state in the (/i, C)-phase diagram may appear rather 
counterintuitive at first sight. Indeed, due to the presence of the periodic potential, it 
costs less energy to form Cooper pairs with a momentum that lies parallel to the periodic 
potential as compared to Cooper pairs with a momentum that lies skewed with respect to 
this potential. To explain why the FFLO-S state can be the ground state of the system, 
we give a qualitative description of the pairing mechanism of the FFLO state in the system 
under consideration. An insightful visualization of this pairing mechanism is to plot the 
Fermi surfaces of up and down particles in momentum space, as shown in Fig. [31 In this 
figure, eight different pictures each show the Fermi surfaces of up and down spins at given 
values of the chemical potentials /x and C,. Around each Fermi surface, a band of size 2A is 
indicated. Only particles that lie within these 'pairing bands' can participate in superfiuid 
pairing. This means that there has to be an overlap between the pairing bands of up and 
down spins for superfiuid pairs to exist. Each picture in Fig. |3] corresponds to a point on 
the phase diagram in Fig. |2l For reasons of clarity, this phase diagram is reproduced in the 
inset of figure |3l where the points corresponding to each of the eight pictures are indicated 
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by small black filled circles. The different pictures are subdivided in three subsets A, B and 
C, which represent distinct physical situations. 

Subset A corresponds to the situation where the Fermi energy of the system lies lower 
than the energy of the top of the lowest Bloch band [Ep < 26). In this case, both the Fermi 



X. 



spherical and the pairing mechanism 



32| . Picture Al shows a spin-balanced 



surfaces of up- and down-particles are approximate 
is roughly the same as in the known 3D situation 
superfluid, where the Fermi surfaces of both particles overlap completely since there is no 
polarization. In picture A2, a nonzero polarization has been introduced by increasing the 
chemical potential ( above a certain critical value. As a consequence, the system has made 
a transition into the FFLO-P state, where the magnitude of the wave vector of the FFLO 
state is such that it partially re-aligns both Fermi surfaces so that a maximal number of 
FFLO-type pairs can be formed. Figure A2 shows the essence of the FFLO state: granting 
a nonzero momentum to all pairs requires energy, but this energy is gained back because 
more bosonic pairs can be formed. The presence of the periodic potential favors the FFLO 
state, since it becomes energetically more favorable to translate the minority Fermi surface 
along the direction of the potential compared to any other direction. This is because in 
this direction the energy dispersion (given by expression (|3])) is flatter than a free particle 
dispersion, which means that it costs less energy to bridge a gap of given size in momentum 
space. In the case of a 3D Fermi gas, the energy dispersion is quadratic in all dimensions, 
so adding the periodic potential lowers the energy for the FFLO state compared to the BCS 
state and the normal state. This explains the significant enlargement of the FFLO region 
in the phase diagram in Fig. |2] compared to the case of a 3D Fermi gas. 

Aside from this energy lowering property for the FFLO state, the presence of the periodic 
potential has additional effects. This becomes more prominent when the density of the 
system is increased and the Fermi energy becomes larger than the energy corresponding to 
the top of the lowest Bloch band {Ep > 26). An example of this situation is shown in subset 
B in Fig. [3l Here the value of fi and hence also the density of the system is larger compared 
to subset A. This can for example be seen in picture Bl, which represents a spin-balanced 
superfluid. In this picture, both Fermi surfaces are no longer spherical because they have 
'hit' the edge of the first Brillouin zone. When polarization is introduced, the system again 
becomes of the FFLO-P type, as shown in pictures B2 and B3. Here it becomes clear why 
the wave vector of the FFLO state increases when the polarization is increased. In order 
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FIG. 3: Visualization of the pairing mechanism of the FFLO state in the presence of a ID periodic 

potential for different values of the chemical potentials. In each picture the Fermi surfaces of up 

and down spins are shown in momentum space. Around each Fermi surface a band of thickness 2A 

is indicated. Only particles that lie in these bands can form superfluid pairs. The corresponding 

points in the (/i, C)-phase diagram are indicated in the inset. In subset A, the Fermi energy lies lower 

than the energy corresponding to the top of the lowest Bloch band. In subset B and C, the Fermi 

energy lies above the energy of the top of the lowest Bloch band. When polarization is introduced, 

the system becomes FFLO-P and the minority Fermi surface is translated by an amount Qz (A2). 

When the density or the polarization increases, the wave vector of FFLO increases likewise in 

order to keep both Fermi surfaces aligned (B2). The extreme case is shown in B3 and C2, where 

Qz = Ql- At even higher density (subset C), the FFLO wave vector becomes skewed with respect 

to the direction of the periodic potential (C3). This costs more energy to translate the Fermi 

surface, but more pairs can be formed because the overlap between pairing bands is increased. 

This lowers the energy enough for the FFLO-S state to become the ground state. 
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to maximize the overlap between the pairing bands of up and down spins, the bulge (in 
the kz direction) in the Fermi surface of the down spins must align with the waist (in the 
kz direction) of the Fermi surface of the up spins. This requires that the Fermi surface of 
the down spins be translated increasingly if the polarization is increased. In picture B3, 
the polarization is large enough so that the FFLO wave vector becomes equal to the wave 
vector of the periodic potential, as indicated by the arrow in this picture. This is the largest 
possible value that the FFLO wave vector can have. Any increase in the magnitude of this 
vector would simply be equivalent to an equal decrease in the magnitude of this vector. 

In both subsets A and B, the FFLO wave vector always lies parallel to the direction of the 
periodic potential. In subset C however, a case is shown where this no longer applies. Here 
the density has increased even more compared to subset B, as can be seen for instance in 
picture CI. Picture C2 shows another FFLO-P case, where Qz = Ql- When the polarization 
is increased compared to this case, the majority Fermi surface will widen and the minority 
Fermi surface will narrow further. This will inevitably result in a decrease in overlap between 
the pairing bands of up- and down-spins, thus raising the energy of the FFLO state compared 
to the normal state. As can be seen in picture C3, the FFLO state can solve this problem by 
forcing its wave vector to lie skewed with respect to the direction of the periodic potential. 
Now it becomes intuitively clear why the FFLO-S state can be the ground state of the 
system. Although it costs more energy to add a perpendicular component to the FFLO 
wave vector (since the component Qz remains equal to the laser wave vector Ql), more 
superfluid pairs can be formed because the overlap between the Fermi surfaces of up and 
down spins is increased, relative to the case where |Q_l| = 0. This is only true when the 
density is large enough so that the Fermi energy lies above the energy of the top of the lowest 
Bloch band [Ep > 26) and when Qz = Ql (subset C in[3]), which explains why the FFLO-S 
state only occurs on that part of the (/i, ^)-phase diagram. One can see the FFLO-S state 
as a final straw at which the FFLO state will grasp to maintain its existence in a highly 
polarized Fermi gas. If the polarization keeps increasing, the cost in energy will eventually 
become too high and the normal state will be energetically more favorable. 
Experimentally, the FFLO-S state can provide an additional signature for observing FFLO 
superfluidity, in comparison with the FFLO-P state. When attempting to observe the FFLO- 
P state, it may be hard to distinguish this state from the ID periodic potential. This is 
because both share the same Bragg peaks in momentum space when the FFLO state is 
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optimally enhanced [33|. The advantage of FFLO-S for experimental detection is that it has 
an additional momentum component that lies perpendicular to the momentum component 
of the periodic potential. Observing this extra momentum component can provide direct 
evidence of polarized superfluidity of the FFLO type. 

IV. CONCLUSION 

In this paper, we have studied the effect of a ID periodic potential on the pairing mech- 
anism of the FFLO state in a 3D Fermi gas. Starting from the partition sum of the system, 
the free energy was derived within the saddle-point approximation. By choosing a suitable 
saddle point, we included the FFLO state in our description, were the FFLO wave vector 
was allowed to lie in an arbitrary direction. Minimizing the free energy with respect to the 
band gap and the FFLO wave vector allowed to study the different competing ground states 
of the system. This subsequently led to the phase diagram as a function of the total and the 
imbalance chemical potential. From this phase diagram we have found that, surprisingly, 
the wave vector of the FFLO state can lie skewed (FFLO-S) with respect to the direction 
along which the periodic potential lies. Furthermore, we showed that this FFLO-S state 
only occurs when the Fermi energy of the system lies above the energy corresponding to 
the top of the lowest Bloch band. To gain a deeper understanding of this counterintuitive 
phenomenon, we explored the pairing mechanism of the FFLO state by studying the Fermi 
surfaces of the spin up and spin down fermions in momentum space. This visualization 
showed that the FFLO-S state can be energetically favorable with respect to the FFLO-P 
state, because tilting the FFLO wave vector can result in an increase in the overlap between 
the pairing bands of up and down spins. We argued that, experimentally, the FFLO-S state 
can provide an additional signature for FFLO superfluidity in a 3D Fermi gas. This is be- 
cause this state has an additional perpendicular momentum component, which allows it to 
be more easily distinguished from the ID periodic potential than the FFLO-P state, since 
the latter shares the same Bragg peaks with the potential when the FFLO state is optimally 
enhanced. 
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Appendix A: Analytic treatment of the free energy 

To the best of our knowledge, it is impossible to calculate all integrals in expression ([8]) 
for the free energy analytically. It is however possible to solve the radial part exactly. To 
accomplish this, we introduce a cutoff K^. in the integral over k^ and divide expression ([H]) 

into three parts: 

1 IA|2 

(Al) 





^l-sp — 


1 

(27r) 


3^1 


+ h 


+ h) + 


A' 
Stt 


where the three different 


parts are 


given 


by 









' ^1 = C dO 1% dK {i:- +... + /,t^ + C') dk^ k^ max (E,,q, |C.,q|) 
h = C dO IX dk, j- dk^ k^ (Ek,Q - ek,Q - ,|,^^,^,'!l(^)J| ) , (A2) 

^ /3 ^ Cdojx dk.j-^'dk. k. (-a,Q - ,|,,,4,'l(^^j| " 



and the shorthand notations iQ - (TTT!) are used. In the first part Ji, the different values 
fci, ..., k^ represent the roots of the equation -Ek,Q = iCk.ql, solved with respect to k^. These 
roots are solutions of a fourth order equation and are rather cumbersome to treat analytically, 
so it is best to calculate these values numerically. The main purpose of the division in (IA2p 
is to separate the logarithmic divergence in the integrand in I3 (when k^ = 0,k± = 0) so 
that this part can be treated analytically. The terms in I2 have to stay together because 
they all diverge for fc_L — )■ 00, and they only cancel out when combined. We now proceed to 
treat each part analytically to the furthest possible extent. 

To calculate the integral over k±_ in Ji analytically we consider two distinct cases. If 
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-E'k.Q > ICk.ql the result is given by 



h 



2n 



(19 



'Af+2N 



dh[Yl 



i=l 



X < ^ log 



[Kl, + A {K)] + J [Kl, + A (fc,)] + A2 






(A3) 



where the set of values Ki = {0,ki, ...,kN, Kc) consists of all roots of equation -Ek.Q = |Ck,Q|, 
complemented with zero and the cutoff K^. In the other case when -Ek,Q < ICk.ql the result 
is 



h 



2-K 



dO 



'N+2\ 



dK\Y, 

+ ((1 - 5Qxcos(e),o) {© [Q± COS [9]] - [-Q^ COS {9)]}) 
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(e \K, - knp] - e [kr^p - K,+^]) 
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+9 [/c„p - i^i] e [Ki+i - knp 
Q± COS (9) 



[^i+l '^^np + -f^i j 



{Khi - ^kl + Kf) 



(A4) 



where the following shorthand notations were used 

S(y=C-«=in(^)=in(^| 



/^ 



Q? 



fcn 



-5 (fc.) 



(A5) 
(A6) 

(A7) 



"^ Qx cos (9) 

In both expressions flASP and flA4p the set of values i^j = (0, /ci, ..., fciv, .K'c) consists of all 
the roots of equation -Ek,Q = |Ck,Q|j complemented with zero and with the cutoff Kc- 



16 



Expression I2 can also be partly calculated analytically. Here we simply state the result 



/ 



TT 



Ql 



Ql 



dkr 



A^K'c + Aih)] 



[Kl + A {K)] + ^[Kl + A {K)f + A2 
2{Kl + e{K)) 
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-A^ log 
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with 



nk 



£ik,) = S 



1 — cos 



Ql 



(A9) 



The last part Is contains a logarithmic divergence in the integrand, but fortunately we found 
that this expression can be calculated analytically: 

Ql 



h = --kK'Ql <K'+2 



+ QL\A\n\o^ 



6- /i- 



6 



{Kl + 5) + y^Kp^2lQ5 ' 
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